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ABSTRACT 


Large amplitude oacillationa of strings under forcing induced 
by transverse loading and the axial movement of supports, exhibit 
several features near their resonance frequencies. The response 
characteristics of such multi-parameter systems undergo drastic 
changes as the parameters are varied in the neighborhood of a 
certain set of their values. For these values of parameters, the 
governing equation of motion is singular. Under such conditions 
more than one set of response characteristics coincide. The system 
follows the one which is stable. This information is represented 
graphically as bifurcation diagram. Bifurcation diagrams are plots 
between the response variables and the parameters. 

In this work, an attempt has been made to study the 
stationary solutions for the system when transverse forcing is near 
a resonance frequency, and movement of supports is done in the 
neighborhood of twice that frequency. The bifurcation diagrams for 
the conditions for which the governing equations are not singular 
have been determined by direct inversion of the equations. This 
procedure fails as the parameters approach the closed set of values 
for which the evolution equation is singular. For such eases, the 



singularities of the problem are classified on the basis of their 
bifurcation characteristics. This property helps in representing 
the possible bifurcations around these critical values in minimal 
terms and also is of adequate use in the study of multi-parameter 


systems . 


Chapter 1 


INTRODUCTION 


The analysis o£ non-linear evolutionary problems is quite 
critical for physical systems, mechanical or otherwise. Non-linear 
equations are frequent in mechanical systems whenever an attempt is 
made to model the systems more precisely. Even in the domain of 
linear physical relations, like stress-strain relationship, 
geometrical and other system specifications coupled with Newton’s 
laws of motion lead to a set of non-linear evolution equations. 

As typical of even simple non-linear equations, non-linear 
models may predict one or more steady state solutions. These 
solutions correspond to the system configurations, after lapse of 
reasonably large time interval. If certain parameters are to be 
varied in the equations of motion, there may be even more 
complexities such as new emerging steady state solutions. In 
physical sense it may become possible for the system to change its 
response characteristics. The physical systems, when confront such 
situation, follow the solution which shows temporal stability. In 
other words, those solutions are prefered to which all the 
transient solutions in the neighborhood converge with passage of 
time. Hence, study of such situations, when many new steady 
solutions emerge from one such solution on variation of parameters, 
are critical to decide the actual response of the system. 

In the present work an attempt has been made to analyse the 
equations of motion of metallic string under transverse and axial 
loading. The response of a string to such forcings is the simplest 
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example of continuous systems exhibiting characteristics not 
explainable by the linear theory. 


11 EXPERlhtEm'AL RESULTS FOR STRING UNDER TRANSVERSE LOADING 

It was first reported by Nur-r-isan. CiQ48> that when a metallic 
string is placed in a magnetic field and alternating current is 
passed through it, even though excitation is always planar 
(perpendicular to the magnetic field) conditions exist under which 
the planar motion becomes unstable and wire begins to whirl. 

This phenomenon is a consequence of an internal resonance due 
to presence of non-linear coupling terms in the equations of 
motion. Due to symmetry of the system, the frequency of motion in 
the plane perpendicular to the plane of excitation is same as that 
in the plane of excitation. Thus, the two motions are strongly 
coupl ed . 

The behavior of strings under different sets of conditions 
has been an area of much experimental interest. L€f& CiQSl^t 
Ramalvc'tshxia. and Wiirthy GiQGSy» Ke^'at and Mctthj&ws CiSUD and 
Ell&x- Ci&lSD conducted experiments under different sets of 
conditions and reported phenomena of jump and whirling. 

General observation of the experiments conducted on the 
strings under planar excitation is that, response of string is 
planar provided the response amplitude is smaller than a critical 
value. Above this value the planar motion ceases to be stable and 
gives way to whirling motion. 

As the frequency is varied from a value far below the natural 
frequency to a value far above, the observed response of the string 
is at first confined to the plane of driving force, but as soon as 
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amplitude of response exceeds a certain value, the planar motion 
becomes unstable and the string acquires an out-o£-plane motion. 

As the frequency is increased further, the amplitudes of both 
in-plane and out-of-plane motions increase until a critical 
frequency is reached at which the out-of-plane motion disappears 
and the amplitude of in-plane motion simultaneously becomes 
smaller . 

On the other hand, if the frequency is decreased from a value 
far above resonance frequency, the response is first confined to 
the plane of excitation until a critical frequency is reached at 
which the response amplitude becomes much larger and response 
string is no longer confined to the plane of excitation. As the 
frequency of excitation is decreased further the response amplitude 
decreases, and when it goes below a certain critical value the 
out-of-plane component vanishes. 

The frequency at which the upward jump takes place is less 
than that at which the downward jump takes place. This indicates 
the hardening of string due to non-linearities. 

12 ft£V2£V OF MATHEMATICAL MODELLIN6 EFFORTS 

Attempts have been made by the researchers to model the 
motion of the string incorporating such behaviors. Anarui CiQ6S:> 
studied the problem incorporating viscous energy dissipation. His 
analysis showed that damping is critical for any such model. Only 
when damping is taken into account, it is possible to predict the 
frequency at which resonance or jump occurs. 

Incorporating viscous damping in the equation of motion 
provides an important conclusion that resonance frequency of a 
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non-linear string is a function of driving force. Resonance 
frequency increases with the driving force. Increasing 
non-linearity of the string leads to increase in the resonance 
frequency and decrease in the response amplitude. 

The effort in this work is to find and analyse the more 
complex phenomena that may occur as a result of combination of 
transverse forcing , and parameteric excitation due to movement of 
supports, in metallic strings. 


1.3 SYNOPSIS OF MATHEMATICAL PROCEEURES AND TECHNIQUES FOR 
ANALYSIS OF CONTI NUOJS SYSTEMS 


Study of large amplitude vibration of continuous systems 
provides clues to their response characteristics, which are not 
explained by the linear theory. Strings, beams and tubes are in a 
way similar as they can be catagorised as one dimensional 
continuous systems. They show behavior quite similar to those 
discussed in the previous section. 

Plates under excitation also show response characteristics 
which are not explained by the linear models. Flutter is one such 
phenomenon in which two dimensional continuous objects interact 
with fluid flow. 

Critical nature of response of non-linear systems pose 
significant problem in design, especially with the design of 
dynamic systems. Stability analysis is an integral part of any such 
effort . 

For systems governed by non-linear evolution equations, exact 
solutions are not generally obtainable, hence recourse is to 
approximate analysis. The methods used for analysis generally 
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consist of purely analytical techniques, purely numerical 
techniques, and numerical perturbation techniques. The purely 
analytical techniques are only applicable to the systems with 
simple geometries and boundary conditions. Purely numerical 
techniques such as finite difference analysis or finite element 
analysis require large computation capacity and time, and do not 
bear much insight about general nature of response of the systems. 

Analysis in this work consists of numerical perturbation 
techniques and attempts have been made to analyse the problem 
taking some ideas from BIFURCATION THEORY and SINGULARITY 
THEORY. 

As the non-linear problems do not allow exact solutions 
recourse is often, to try to study the solution in which some 
characteristics remain invariant with time. Physically, this 
corresponds to analysis of the general behavior of the systems over 
a large time interval. These characteristics may be coordinate 
invariance in some linearly transformed domain C stationary 
solutions , periodicity or a combination of both. BIFURCATION 
THEORY is a theory of such equilibrium solutions. The theory 
discusses the methodology by which such solutions may be 
identified, or the conditions required for existence and stability 
of such solutions may be prescribed. 

In non-linear evolution problems of many variables and 
parameters the solutions in themselves contain little information 
about general nature of the system, it is the interaction between 
the solutions of different types that provide clues to the general 
behavior. Such interactions are possible at the singular points 
i.s., the points, where the solutions intersect or linear terms in 
the determining equations vanish and the solutions around these are 
determined by non-linear terms. 

For this reason it makes sense to classify the solutions by 
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their organisational properties such as type of bifurcation that 
the solution is capable of displaying with variation of parameters. 
Based on such a classification, implications of variation of a 
parameter, or a combination, on the nature of the response may be 
predicted. This process of classification comes in domain of 
SINGULARITY THEORY. 

14- GENERAL NATURE OF CONTI NUOS SYSTEMS 


Continuous systems such as strings, beams and plates are 
known to show certain phenomena, as whirling in strings, beams and 
tubes, and flutter in case of plates. These phenomena are not 
explained by linearised equation of motion. Non~linear behavior of 
response in such cases can be attributed to large amplitude 
vibration in the system. 


For small amplitude oscillations, response can be adequately 
described by linear equations and boundary conditions. At higher 
amplitude the non-linear behavior becomes more prominent. Uhirling 
and Jump occur only at amplitude more than a certain value. Motion 
evolution equations modelled on large amplitude deformation in such 
systems have higher order non-linear coupling terms which become 
significant at large displacement. 


The source of such non-linearity in 
geometry of the system. At higher amplitude 
becomes a non-linear phenomena and leads 
relationship between strain and displacement. 


the string is the 
stretching of wire 
to a non-linear 


To facilitate the analysis the non-linear couplings are 
considered to be a perturbation from the linearised system. The 
solutions of the linear equations are used as spatial variation. 
Using orthogonality of the mode shapes, a non-linear coupled 
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second order differential equation is obtained. This describes the 
temporal behavior of the system. 


The latter equations are reduced to a set of 
differential equation using a perturbation technique such 
method of multiple scales. The equations so obtained 
proved to be an approximate description of the actual 
process . 


ordinary 
as the 
can be 
physical 


1.5 THE METHOD OF MULTIPLE SCALES 


The technique of multiple scales [2] consists of viewing the 
solution which is a function of one variable C time, t ), as a 

function of many variables, namely ( t, «t , «*t, ,«^t ), where « 

is a small positive number. In this method the expansion 
representing the response is considered to be a function of 
multiple independent variables, or scales, instead of a single 
variable, the process starts with introduction of new independent 
variables given by, 

T = «’'t for n = 0, 1, 2, (1.01) 

Tk 

Thus, the derivatives with respect to t, becomes expansions 
in terms of the partial derivatives with respect to T^. 


d dT d ^ dT d ^ 

O * 


D 


o 


+ «D 

4 . 


+ 



= d“ + 2«D D 

o a 4. 


c“(d“ + 2D D ) + 

' 4 0 4 


( 1 . 02 ) 


The solutions are assummed to be expansions of form. 


:(t : 4S) =«x^( Tq, T^, Tg, ) +«*X2^^0’ ^1’ ^2 ^ Cl. 03) 
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In this expansion the number of independent time scales 
needed depends on the order to which the expansion is carried out. 
If the expansion is carried out to order of then T and T 

O A 

are needed. For example, consider the system governed by equation 
of the type 


. . M n 

X + S a X =0 (1.04) 

n=4 ri 


3 

The expansion carried out to 0(« ) requires the first three 
terms of the expansion for and x(t; «) . Substituting the 

series in the equation and putting the coefficients of «, to 

3 

be zero so as to satisfy the equation upto order « , we get 


0(«) ; 

n2 ^ 2 

D X + X = 

O A O A 

2, 

2 

2 

0(« ); 

D X + X 

St OS 

3, 

2 

2 

0(« ) : 

D X + <.> X 

os OS 


= -2D D X X 

O A A a A 


(1.05) 


2 3 

= -2D Dx-D X-2D Dx-2oi xx-oi x 

OAB AA OBA BABBA 


The next step is to substitute the solution of equation 
(1.05. £) in equation <ri.05.ti) and to look for secular or resonant 
terms in the resulting equation and equate them to zero. This 
process makes x^ to always be of order abx^ as, the solution of the 
resulting equation with all the resonant terms zero, will always be 
bounded. 


This reduces the need to study the resulting equation as it 
has no bearing on the general non-linear behavior of the solution. 
Solving this equation would obviously help in increasing the 
accuracy of the prediction. 

The coefficients of the resonant terms when equated to zero 
give an evolution equation for x^ in terms of large time scale T^ . 
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The evolution equations so obtained govern the variable of order « 
in large time scale. These equations hold clues for non-linear 
behavior which occurs in the system. 


t6 STATIONARY SOLUTIONS OF EyOLUTION ECUATIONS 


The evolution equations in the large time scale for response 
of order e can be represented in following form. 

— =FCx,X) ( 1 . 06 ) 

For steady state solutions of coordinate invariant type, 
these equation symbolise a mapping from one finite dimensional 
space to another. 


F : « X A ► V 

F^: « > V 


(1.07) 


The stationaS^ solution of the equation (1.06), for any fixed 
value of X , can now be observed to be the set of intersection of 
the graph of the map (1.07.i£.> with V = 0 in 5^ x W . The 
parameterisat ion of the set by the set of parameters ( X ) is 
represented graphically as bifurcation diagram. 


As the aim la to construct and study such diagram to predict 
the response characteristics of the system as they change with the 
variation of the parameters, it Is of interest to study the 
solution set in ( x A ) about a point (0,0). 

The Taylor aeries expansion for any such mapping with 
8 * i ~ (0 < «5 << 1) corresponds to 

V = D F(0 ,X)x + D®F(0 ,X) + (1.08) 

M X 4 
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For the c&ses when D^F(0,Q) is invertible, introducing s 
coordinate transformation. 

U = D 

Taking * <5, the series is reduced to following form. 

U = X + G(x,X) , (1.10) 

Where |j D^G(x,X)| ~ <5 
As for our purpose we require U = 0, hence 

+ G(x^,X) = 0 (1.11) 

o o 

This equation can be solved for x in terms of X by a simple 
procedure and solution in the neighborhood of (0,0) can be 

approximated by a form containing X only. 

"^0= " G(x^,X) 

X® = - G(0 .X) 
o 

X* = - G(x^,X) (1.12) 

o o 

X*' = - G(x’'~*,X) sa X (X) 

This process converges, ie , valid if x and X are 

sufficiently small. This procedure leads to an approximation of the 
bifurcation diagram about a point (0 ,0). It is due to 

invertibillty of D^F(0,0) that it has been possible to obtain 
explicit representation of response as a function of parameters. 
Now, consider another change of variables in the equation (1.10), 

Y = x - x^(X) (1.13) 

Taking | Y | << 1 , following form is obtained. 

U = Y (1.14) 


F(0,0) 


(1.09) 
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The graph o£ the equation (1.08) and the map (1.14) represent 
the same zero set. 

As an observation, it can be contended that any expansion 
having invertible Jacobian can be mapped on such a form (1.14). As 
the set of invertible matrices is open, for any small change in the 
parameters, the Jacobian remains invertible. Thus, in the 
neighborhood of X = 0, the solution is directly paramet erised by 
the parameter set, X, and does not show any bifurcation. 

This procedure can be used to study the bifurcation behavior 
for a general map. The representation (1.14) has been obtained by 
linear transformation on the relevent coordinate system. Hence, 
bifurcation properties of such representation is same as those 
shown by the maps it represents. 

The non-linear nature of the problem allows invertible 
transformations in the coordinates. This property is used to 
reduce the problem to minimal dimensions. Analysis based on such 
theory helps amply in plotting the bifurcation diagram as due to 
non-lnvertibility of Jacioblan in the neighborhood of singular 
points, simpler computational procedures fail. 
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Chapter 2 


FORMULATION 


The equations governing the motion of string are determined 

by Newton’s laws and, geometric and material properties of the 

system. Formulation of the equations of motion is based on the 

same presented by Nayfeh and Nook [2]. To obtain the equation of 

motion , we consider the equilibrium of forces on the differential 

element PQ (Fig. 1). At any instant of time , consider a vector 

field (u, V, v) which is a function of time and spatial variables. 

The displacements are measured from the dynamically undeformed 

position, and x-axis is setup along the string. The string in its 

undeformed position is assummed to be taut and has tension of 

value hi . 

o 
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The length of the differential undeformed segment of the 
string PQ is taken as dx. Let the coordinates of such a segment be. 


r = X I 
p 


r = (X + dx) i 


( 2 . 01 ) 


Then, under a displacement field (u, v, w) , the displaced 
coordinates of the points P and Q are. 


,= (X + u) f vj -f W fc 


( 2 . 02 ) 


r^,= (x + u + dx + u,«dx) * + (v + v^dx) + (w + w^dx) is. 


The relative displacement vector r ^ is given by. 


r . , = (1 + u)dxf + vdxj‘ + wdxM 

p ’iqi * K K K 


( 2 . 03 ) 


The length of deformed segment ( PQ ) at some 
(t, X, y, z) is given by. 


point 




(1 + u +v* +w * dx = ds 

X K K 


( 2 . 04 ) 


Thus, the unit tangent vector to the deformed string segment is. 


(l+u)i+v j' + 




( 2 . 05 ) 


The instantaneous value of tension in the deformed segment of 
string is then given by. 


ds - dx 


N = N + EA(- 
o 


( 2 . 06 ) 



Using the above information equation of the 

differential segment can be formulated as, 

|^[ N ^ ] = pA C ) - F (2.07) 

where, 

p is mass per unit length 
A is area of cross-section of string 
F is external force per unit length acting on the 
differential element. 

Using equation (2.05) and (2.06) , equation (2.07) can be written 
as , 




[ 


EA)^ + EA)((1 * xxj t * J + 


M)j 


pA( i + v^^ j' + w^^ k) - F 


(2.08) 


In x-direction the force balance is represented by, 


pAu - EAu = (N - EA )^ 

It HH O VX 


r (1 +u ) 


[(1 +u ) +v +W^ 3 


+F 




r ^ ''x - - 1 

fd + ) * 1 

L (1 + u ) ■* 


(2.09) 


+ F 


At this instant another assumption is made, that amplitudes 
are small compared to the possible wave-length of the modes; hence 
the approximation. 


2 

V 

X 


w* « (1 +u )* 

X X 


( 2 . 10 ) 
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I 


This approxifflstlon reduces the equations to following forms : 




= - CEA - N ) 

a a' 


[-[ 


Cl -2u )(v + 

H X 


<^]] 


( 2 . 11 ) 


For the governing equation in a transverse direction (y), 
pAV -EAv = (N - EA)f — f a 1 

^ ti KK o ■'ax I r ^ a, — I 

L[(l+u)+v +w] aJ 


r ^ 

I v - , 

|_ XX 


(N^- EA){v - -^(v (u 


a la 1 a 1 

u + — V +-rw))i 
X 2 X 2 X "J 


pAv N V 

tx O XX 


= (EA - N )^(ev ) + F 


<i'''Sx^ X y 


(2.12) 


where , 


e = u -u* +-i v“ ^ 

X X 4 X Z X 


Similarly, for z-direction momentum balance , we have 


pAw^ -N w = (EA - N )|- (ew^) + F 

it O XX Cr &X. X ■ 


(2.13) 


In the formulation , 


■ai 


= [(1 +u ) +v +w r • 

** ' X X X 


. a 1 a 1 a 

= l-u +U - — 'V - — W 
X X a X ax 


= 1-6 
ds-dx 


e = ( 


ds 


(2.14) 


Thus, 6 denotes strain in the string . By symmetry of the problem e 
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can be taken to be a periodic function of time. Putting, 


c; = E/^ 


C = N /pA 

B O 


the set of governing equations becomes 


2 2 2 

= (C. - C.) [ (‘ -«^)Cv; 


+ F 


V - C V 

t.t S XK 


cc* - e-,t-\ T (u V + i V* . i .1 

‘ <> .kJ 


tt 


C*w = (C“ - c“)J_r w cu -u* + i V* . 1 >1 

axx d adx[KHx ax axj 


(2.15) 


+ F (2.16) 
y 

+ F 

z 


For metallic strings it is known that 
C®/C* = 0(400 - 1000) 


(2.17) 


Ue now non-dimeniional ise the set of equations (2,16) and 

assume the forcing to be of order one. Then, the amplitudes in the 

transverse directions are of order one. From second equation, the 

amplitude of longitudinal direction oscillation is of the order 

C*/C*- This result when interpreted for the first equation means 

2 

that order of magnitude of the inertial term is 1/C^ times less 
than that of other terms in the equation ( 2.16.i ). Hence, due to 
these considerations, the longitudinal inertia term may be safely 
omitted . 


- C u 

d KX 


= - 




S 

Sx. 


f (J- -u )(v* +w“)l 

L “ ** “ J 


(2.18) 


Here we assume dynamic forcing F^ = 0 
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Integrating the equation ( 2.18 ), 

2 2 

-u + «(t) = - u (v“+ w“) C2.19) 

X -S— K K K 

This equation makes it clear that u is of same order as v and w . 

H K 

fi S 

Hence, to the order of (v , w ) , 

u = «Ct) - — (v® +w*) (2.20) 

K ^ K K 

An important implication of this result is that, 

eCt) = oi(t) (2.21) 

The equation ( 2.19 ) is integrated again, considering following 
boundary conditions on u. 


x=0 


= 0 ; 


^ix=/ ' 


( 2 . 22 ) 


where p(t) is taken to be a periodic function. The integration 
yields , 


X 

u = x6(t)-i|Cv* + vf)dx 
2 Jo ** 


Hence, 


e(t) = p(t) 


i 


(2.23) 


(2.24) 


Adding viscous damping terms 
the transverse direction equations. 


and eleminating e(t) from 
the system of equation reduce 
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to. 


ti 


+ 2^v. 


-C“v 

ft M 


w +2/JW -C V 

it 1 ft XX 


C*PV + 

4 XX - ■ 

X 


cfPw + 

A KX 


C“w 

A XX 



+ W®) 

X 


+ W ) 

X 


dx + F 
y 


dx + F 

■ 


(2.25) 


(2.26) 


2.1 LJWE^ft WOOES APPROXIMATION 


As V and v are continuous and bounded for all x, they can be 
approximated by following series expansion using linear modes 
solution as spatial variation. 


1 

— <sc 

V = (t) Sin(n»Tx//) 

Ti 

ri«A 

(2.27) 

4 

— oc 

w = ct) Sin(nnx/-f) 

^ ri 

ri-* 


± 

— OC . 

'“’f EC, 

A 


Sin(n*tx/-^) 


± 

— oc BBSS 
+«*■« EC C « Sin(nf*x/^) 

F» ft 

ri«4 


A 

— OC , 

n*A 


SinCnrr x/ip ) 


A 

•*- CA ft ft ft 

= - c“ P «“/ E" " C^sin(nnx//) 

r»-A 


3 

— ftOC ft ft ft 

(f E C^»' Sin(ri«x/-^) 

r»M A 


IJ- 


OC 

E 


(f mn/'/ 

Trt 


Cos(mnx/^)) + ( 


oc 

E 

ro*o 


n SBTI /X 

fW 


Cos(Hnx/Jr) ) 


+ e(x,t) (2.28) 
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Since , 


Aa, 


I 


C/C 


oc r ft ft s ft ft 


( E C mn// CosCmfTx/^)) dx 


O m»»0 




vx n /i Cos Cwix/<f)dx 

m 


L 


CosCmnx//) CosCnnx/JT) = •»• *5 

z mn 


(2.29) 


(2.30) 


The integral ( 2.29 ) is equal to. 


Q£ 


= EC 


ft ft 
ft m f7 


m 2 ^ 


(2.31) 


Thus, using orthogonality of linear inodes, the set of equation 
reduces to the following form. 




C + -2-^ C“c + 2 #J K 


ft ft ri 


-(f 


ft ' fS 


c“«“^c — ^ E ( «“cc*+ r)“ )) + g(t) 

4. rk - ^ mm 

4 « m 


(2.32) 


Nov, define some new variables 6> , F , p denote 


2 2^2 4.2 

„• = £_!!_=. , r = ^ 

" i* 4i* 


IT 


*0* 


; «p = P J.* ; = fj ; 

4r 


(2.33) 


Introduction of these variables reduces equations to the form 
given below : 


C + « C = 

Ti rt Ti 


^ n-pc^ *r-nc„i: C”(C' *»*>)] 


+ & (t) 
®h ' 


(2.34) 


YJ + <■> T) 


= - «r2/j + n“p Y) + rnn E Cs? ( C* + T9* ))1 + f Ct) 

I ri n r* m wi I n 
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22 PERTURB ATI ON ANALYSIS IN THE NEIGHBORHOOD OF A RESONANCE 

FREQUENCY 


Next step consists of analysing these equations in the 
neighborhood of resonance. As apparent from the equations, 
non-linear behavior becomes more prominent in its vicinty. Let, 


g (t) = 2 ts k Cos Ot 

Tt r» 

f (t) = 0 

p(t) = 2 S Cos (2ft t + y) 

p 

ft = 

r* 4 

ft 

p = 

r» a 


(2.35) 


At this stage the method of multiple scales is applied. 
Consistent with the method, the solution is assumed to be an 
expansion in terms of the small number 


C (t ;« ) = C (T ,T ) + 4 bC (T ,T ) + 

n ’ ' ^ no 'o’* rsi'o’* 


y) (t ;e ) =>)^CT ,T) +«n.(T .T) + 

ri rvO O 4 ri4 O 4 


For order 


D*C + C =0 

O r»0 ri nO 


Y) + n =0 

O no n no 


(2.36) 


(2.37) 


The first order 


C ^(T . T ) = 

no O A 


approximate solution is given as, 

A (T )exp( iu T ) + A (T )exp(-i6> T ) 

ts 4, irsO n4 rs 


n (T , 

no O 




B ( T )exp( i6> T^) 

rt 4 r% O 


+ B^(T^)exp(-i4^^T^) 


(2.38) 
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For order 


D C f = 

O r»4. n ri4 




+ 2k Cos (ft T ) 

TV O 


(2.39) 


+<^*T7 = -2D D 

O n* n r»4. O 4 


"no-2f'D„n„-pn\^-n*r„_^j^ £ ■"’ (cl. "n!. )l 


mO fwo 


Putting the solution ( 2.38 ) in the equations (2.39 ), as 

required by the method, following form is obtained, 

■>X *“X= -2'>«D a e-‘“„To -20 (A ^".^0 *A ) 

” o rt IS 

1% rs 

-rn* ( )f E n? ( (J^ b?c, -fA e’^^ n?^ f 

L 

■•• ( B„, e~’'“m )“ ) j + k^ ( ) (2.40) 


To eliminate secular terms from the expansion ( 2.40 ) the 

coefficients of e rv o and e ' n^o have to bee zero, 

-2i.i.i (A +pA ) -n “« e '■'^A e% 1 -2rn“Ar in“(A A +B B ) 

" " ri r, ' f« wi mm 

Tn 

-rn-A^CA^B') = 0 (2.41) 

where ’ denotes derivative with respect to T . 

For similar reason. 

Complex conjugate of the above equation = 0 
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For the other transverse direction, 


- 2 l 6 > (b’ +ph ) - 


2 rn*BX; m* (A A +B B ) 

ri*" mm mm 

m 


-rn'‘B (A*+B*) = 0 

ri ri ri 


(2.42) 


Complex conjugate o£ the above equation = 0 


Let , 


A = oi C 

ri r\ f\ 


B = « D 

ri « r 


(a is a scale factor 

Tt 

to be determined later) 

(2.43) 


•2i<a « (C +juC ) -n** a e’‘^c -2 rn*« BT oi*in“(C C +D D ) 

'Ti rs fs f\ rsnn fl ri m mm mm 


-rnVc . 


t\ n ri ri 


(2.44) 


Let the scale factors be. 


= ctr 2 

n 

h = k /(2oi ) 
rs n 

^ = n®* ©'■'^ / 2<J 

ri n 


= Pn /6> 

n 


(2.45) 


Introduction of these scale factors lead to the following set of 
equations , 


-v(C +^C )-^fr C C E «®in (C C +D D )-C ((? )+k e^n^* =0 

nT% rHri""m mmmm ri nn^T\ 

xn 

(2-46) 

“^(D +pD )- 4 i D ©‘■^*'^*-f9 D E «“in (C'C +D 15 )-‘D ((? +lf ) =0 

^ ^ rsry nn m mmmm inri’i^ 


m 


22 



2.3 EFFECT OF DAMPING 


For a mode ( k ) not excited either axially or tranaversally , 

0 

(2.47) 

0 

Then, multiplying equation (2.47.i ) by C 

multiplying complex conjugate of equation ( 2.47.i ) by , 

multiplying equation ( 2.47.£i ) by . 

multiplying complex conjugate of equation (2.47.tt) by -D 

k , 

and, adding the resulting equations, we have, 

■"[ (|Cj*)’ + 2iu|Cj* +(|D^|*)’ = 0 (2.48) 

For non-zero p , this condition can only be satified by, 

|C^| = |D^| = 0 (2.49) 

This means that for mode (k) having no excitation the 
response is always zero. If excitation is restricted to a single 
mode (mode n), then due to damping in the system other modes are 
Indentically zero. Due to this reason analysis can be restricted to 
the mode having its linear natural frequency near the frequency of 
forcing. 


l(C +^C )-^ C e^^2^4-f3 C Ea^m (C C +D D )-C C(f +lf ) 

n n nri rin“m mmmyn 

m 


-(d’+hD )-i^ T3 D c «** if (CX +DD )-t) (0 +ff ) 

^ n ri nri nri^fw mTWwwi r\ r% r% 

m 
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EVOLUTION EQUATIONS FOR THE FORCED MODE 


let , 

C 


y. e 


io- T 

1 4 , 


lo* T 

D = y e 4 i. ; 

ri a 


— T T 

C=ye44; D=ye44; 

n a ri 4 

&= a - 2<y ; 

a 4 

f L dT 

+(p c^i^oy^^e 1 +2 '■y^Cy^yg+y^y^) 

» --LdT 

+C^i -^<y^)yj^-C^Oy^e 4 +2 ^ yj, C y^ y^ +yj, y^ ) +*• 

» 1 . 6T 

y^ +(jU +i.<y^ )y^- C<^u)y^e 4 -2 *■ yg CV^ +ya y^ ) 

> -- 

y^ + Ca* ) y^ - (^'•)yge *■ +2 i. y^Cy^y^^+y^y^ ) 


C2.50) 

y* ^ + i k = 0 

a 4 a 

y Cy* +y* ) - i k = o 

^i. a * 

- y^Cy®+y®) = o 

+«■ y (y* +y* ) = o 

a i a 

(2.51) 


In this work, attempt has been made to analyse the system 
with ( O = 0), t.e. , autonomus case. 

As , 

y = y and y=y (2.52) 

a 'a -* 


taking. 


y = X + vx 
* * a 


y = X + vx 
a a * 


X 


I X 
a 


(2.53) 


y. = 


X 

a 


V X 

■* 


4>= K 


iv 
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following equations are obtained, 


+(P Sin(v'))x^ -(o +? Cos(v))x^ +x^ (3x® + 3x®+x® + 3x* ) +2XjjX^x^=0 

X +(fJ +? Sin(v'))x -(-<y Cos(v'))x -x (3x* + 3x* + 3x*+x* )-2x x x +k=0 
a a 4A4.fts<4a4ft 

x^ +Cfj -? Sin(v'))x^ -(<!>■ +? CosCyJlx^ +x^ (x® + 3x“ + 3x* + 3x* ) +2x^x^x^=0 
X +(fJ 4-? Sin(v'))x -(-<y+? Cos(v'))x -x (3x*+x* + 3x* + 3x* ) -2x x x =0 

^ >4 Sfi4jEfcSI<44.S-4 

(2.54) 

These equations may be written in the following symbolic form : 

H V 

^ + /(x ; fJ, O', <f , k, v) = 0. (2.55) 

The aim of the further analysis is to study the solutions of, 

If + /^ ( X ) = 0 (2.56) 

o 

as they vary with perturbations in the perameters. 
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Chapter 3 


COMPUTATIONAL PROCEDURES 


The equations determining the motion of the string are given 
in terms of the variables of order e, and large time scale. The 
equations for the case of stationary solutions can be considered as 
a map : 


f ( X, X ) = V (3.01) 

The solutions for the problem are given by intersection of 
this map with v = 0. Bifurcation points are special points in the 
solution domain, at which determinant of the Jacobian vanishes and 
the solution curve in (x, X) in the neighborhood of these points 
are determined by the non-linear terms present in the map. From 
these points more than one branch of solution curves emerge, along 
each of which Jacobian is invertible. Obtaining these points is 
essential to plot the bifurcation diagram, and thereby study the 
steady state response of the system. 

At the bifurcation points, the Jacobian is not invertible, 
which is same as the claim that one or more of the eigenvalues of 
Jacobian matrix are zero. The eigenvalues are the roots of its 
characteristic polynomial. The conditions under which the system of 
equations has a certain number of its eigenvalues zero, are 
specifed by the coefficients of the characteristic polynomial : 

P(«) = +a +a oi* +a « +a (3.02) 

£1 & 4, O 


^6 



Hence, for example, to get a bifurcation point having two of 
ita eigenvalues equal to zero, the set of equation to be solved is 
given as. 


f( X, \ = 0 

a^C X, X ) = 0 (3.03) 

a^C X, X ) = 0 

The set of equation so contructed has m>4 equalities to be 

satisfied. For this reason, the set of variables is augmented by 

allowing m-4 parameters to vary. Thus, given a guess value 

(x^, a , ft ), when (x, «) are allowed to vary ,we try to obtain 
o o o 

(x', ft ") for which the equations are satisfied. 

O O o 

3.i GAUSS NEWTON* S METHOD 


The set of equations ( for example 3.03 ) to be solved can be 
represented symbolically as, 

F( Y ) = 0 where, Y = ( x, OC) 

= F ( \ - F'( Y^ ) AYj^ . 

F'^C Y^ ) F ( Y^^^) = F’*C Y^ ) F( Y^ ) + F’* C :> (. \ ) 

Taking, F( Y, ) = 0, 

k+ A 

F’'*‘( Yj^ ) F( Y^ ) = -F’’C Y^ )F ’( Yj^ ) AYj^ (3.04) 

The procedure runs as follows: for a given value of Y^. 

is obtained from the equation ( 3.04 ) using gauss-elemination 

technique. AY is the direction in which the value of the function 
k 
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approaches zero. 


C3.05) 

The value of i> is chosen so as to minimise the function in 
that direction. For this value of v , is obtained which is 
taken as guess value for the next Iteration. 

The termination condititon for the method may be taken as 
step-size C v ) being less than an error parameter, or the norm of 
the function being less than a given value. 

This procedure constitutes the routine QAUSS_NEWTON. The 
routine is used to find a solution to the given set of equations 
in the programs BIFURCATION- LYAPUNOV and DIRECTION. 


3.2 METHOD FOR CALCULATION OF DERIVATIVE OF AUGMENTED EQUATION 


To calculate the bifurcation points, the function and their 
derivatives are to be evaluated repeatedly. To this end, the set of 
equations ( 3.01 ) are stored as a combination of array and 
records. The smallest entity that is to be dealt with is a term. 
In program BIFURCATION! which calculates the bifurcation point, the 
data structure that is used to store the information about a term 
is called terM_1 • 

It has three fields 

; COEFF stores the coefficient of the term. 

: PWR stores coded information about the exponents of 
all variables and parameters. 

! FUNC contains information about the symbolic 
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functions, such as Sin and Cos, present in 
the term. A boolean code TRN3D is used to 
check for corresponding variables. 

An array of such records, having elements corresponding to 
the terms in the function is used to store functions in the 
program. Such a vector array is called FNX_1 • On similar lines, 
vector functions are defined as a 1-dimensional array of FNV_1 . 
2-dimesional and 3-dimensional arrays of functions are similarly 
def ined . 

Operations of differentiation C DIFFERENTIATE > and 
evaluation ( EVAL 5 neccesary for computation of bifurcation point 
are coded so as, to be performed on such structures. These 
operations make it possible to work with function as basic objects. 


3.3 OPERATIOhlS ON COEFFICIENTS OF CHARACTERISTIC POLYNOMIAL 


For evaluation and differentiation of the coefficients of 
characteristic polynomials, it is neccesary to introduce 
another data structure. It is termed ELEMENT_CD in the program. 
It is an array of another record called COEFF_B, each of which 
denotes the coefficients of the characteristic polynomial 
C PC«) )• These are inputed as functions of entries of an n x n 
matrix ( four at this stage ). COEFF_S ane functions with each term 
an n-tuple at maximum, denoting the combination of entries to be 
multiplied and subsequently added to give values for the 

coefficients. Thus operation of evaluation C EVAI COEFF-JS ) and 

differentiation C EVA! DIFF_COEFF on such structure may be easily 

defined. 
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The coefficients are differentiated in the following way. 


=£ ( where, 1 i i 4) (3.06) 

j l, m I ,»« j 

The symbol A. denotes the entry of the matrix. The first 
L , m 

term under summation in formula ( 3.06 ) is obtainable using the 
code of ELEMENT-jCD- The second term can be computed before hand 
and stored in three dimensional array of reals. 

3.4. CALCULATION OF PROJECTION MATRIX 

The procedures mentioned previos sections give the 
bifurcation points having at least the given number of null space 
dimensions. The next step consists of projecting the map such that 
variables are partitioned to span the null space and its 
compliment. To this end, first the vectors spanning the null space 
of Jacobian matrix at bifurcation point, are computed. The 
compliment vectors are taken to be orthogonal to them. The vectors 
are arranged to form a projection matrix in such a way so that, 

if the null space is m-dimensional , variables x^ , , 

denote the compliment and rest span the null space of the mapping. 

This part of operation is performed in the program PROJ_MAT. 
Given a singular Jacobian matrix by the program BIFURCATION f 
PROJ_MAT computes the required projection matrix. 


3.5 EXPANSION OF THE MAPPING ABOUT THE BIFURCATION POINTS 

The next step in the proposed analysis consists of expanding 
the function about given bifurcation point. This part of the 
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operations are performed in the program LYAPUNOV* The given 
function is a 4-dimensional mapping with B-parameters and 
4-variabl6s. For purpose of expansion no distinction is made 
between the parameters and the variables. The process of expansion 
is performed by translating the domain coordinates. 

The description of the basic entity, term of a function, in 
this program, requires only two fields. 


: c stores coefficient of the term 

; p contains coded information about the powers to 
which each variable is raised. 

Routines are written to perform various operations on 
the functions which are coded using such structures. 


: ADD_TERM 
:MLTP_TERM 

: ADD_FN 
: SCALAR_FN 
: SORT_FN 

: MLTP_FN 

: COMPACT 

: EKPAND_FN 

: CUB9T_rW 


adds a term to a given function, 
multiplies two term and returns the 
resulting term. 

adds a function to a given function, 
multiplies a function by a scalar, 
sorts a given function to resemble a Taylor 
series . 

multiplies two functions and returns the 
resulting one. 

removes any term having coefficient less 
than a given value (EPS) from between any 
two terms. 

translates the origin of the variable 
coordinates by the given position vector, 
substitutes a given function in place of a 
variable in a given function. 


The program LYAPUNOV using above mentioned procedures expands 
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the equations about the bifurcation point : 


V = f( X. X) 

“ f C X + ^ ^ ft ^ 

where, (x^, X^) denote the bifurcation point, hence 

V = e( X, X ) 


( 3 . 07 ) 


This expansion process makes purely numerical valued terms in 
the g(x,X) to be equal to zero, but the Jacobian is same as that 
without expansion. Next step consists of rotation of the variable 
coordinate system using the projection matrix already computed by 
the program PROJ_mat. 

X = P X 

- - X ^ ( 3 . 08 ) 

V = P 8(P X ,X ) 

Let x^ denote the variables spanning the null apace and x 
spanning the compliment of the null space. Similarly, let v be the 
variable spanning the range and v^ the variable spanning the 
compliment of range. 

-^0 ■ 

( 3 . 09 ) 

^x 0) is invertible, by inverse function 

theorem x can be expressed locally in terms of (x^, X) for v = 0 

In other words, 

° *^ft ' ( 3 . 10 ) 
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can be solved for . Multiplying equation ( 3.10 ) with, 

o 


• “>] 




we have. 


Cf €r 


(3.11) 


For I I , I x^l and | X| to be much less than one, 


following procedure helps in finding an expression for x 
in terms of (x^ , X) valid for small values. 

x = G(x , , X) 

« a ' ft ' 

o o ' 


X =0 

0( 
o 


c*. 


(3.12) 


The above procedure converges, as 
I x^ - x^ B ^ I ’ ^ft 


k k Tt. 

X - X 
OL Ol ' 


®x 8 = ' l*/^! » 1^1 ) 1 

ot 


(3.13) 


These operations are done by the routine SOLVE in program 
LYAPUNOV aided by procedures FIRST-APPROX and PROCESS. 

Having obtained solution. 


X = X ( x^ ,X ) 
a m ft 
o o 


(3.14) 


the solution is substituted in the equation ( 3.09.itl ) the map to 
be studied is reduced to. 




(3.15) 
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In these equations , D F(x , 0) is nllpotent. 

3.6 PROCEDURE TO OBTAIN TANGEhlT VECTORS TO THE SOLUTION CURi'ES 
AT THE BIFURCATION POINTS 

Following procedure is used to compute the tangent vectors at 
the bifurcation points. At bifurcation points the main problem is 
to find all the stationairy solution curves around that point. This 
procedure is used to find tangent directions to the curves 
unfolding from that point. This stage of computations is performed 
in the program DIRECTION. 

Let Y denote where ot is the parameter with which the 
bifurcation diagram is to be sketched and F ( Y ) is the reduced 
map restriction to the null space. The directions are 
obtained in the process of solving the set of equations given 
below. 


0 = F ( Y ) 

0 = [ D^F( Y ) ] Z 

0 = I Y I* - C3.16) 

0 = ?* 

0 = 1 Z I* - 1 


Given a vector function the procedure MAKE_FUNC assembles 
the above set of equations. 

If Y denotes n unknowns, the above set of equation are in 
terms of 2n+l unknowns. These equations are solved wih Gauss-Newton 
method presented earlier. The solutions as apparent from 
formulation of the equation set, has all its components zero exempt 
those denoted by Z. It can be claimed that Z represents one of the 
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many possible null directions. 


Different possible directions are obtained for different set 
of guess values for the solution. ? is generally initialised with 
value of order one. Y and 2 , for initial guess are taken to be 
equal and of order one. 

Having obtained all the possible tangent direction a small 
step is taken in each direction, and all the variables in the 
compliment of the kernel are also evaluated using equations 

( 3-14 )- 

This sets of (x , x^ , oi) are used as reliable guesses for 
each of the possible solution emanating from the bifurcation point. 


3.7 COm'INUATJON PROCEEAJRE 


The above procedure gives good guesses for one point on all 
possible solution curves around the given bifurcation point. The 
aim of the procedures in this section is to trace the variables as 
they vary with a parameter cn. These procedures constitute the 
CONTINUATION program. The information available are the bifurcation 
point and a guess point near a solution curve , relative to 
bifurcation point. 

The technique used to construct a solution curve consists of 
first obtaining tangent at a given point in (x, a) on the solution 
curve and moving in that direction by a step-size. The value of 
step-size is governed by various considerations. The most 
important of them is that while moving along the tangent direction, 
the norm of the function does not Increase beyond a preset value. 
Constraints are put directly on the step size, as well. Its 
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maximum value ia set to bring a homogeneity in the distance between 
two solutions on the same curve. Uhen the value of step size goes 
below a certain value the procedure is coded to stop and plotting 
of the next solution starts from the given guess point near the 
bifurcation point. 

Taking step in following fashion, provides a good guess 
values for the next stage of computation. 


Y = (X. «), 

be the point on solution curve known before hand. 

<5Y be the tangent at the point, 
o 

Y = Y + step X <5Y . (3.17) 

4 0 0 

This process is repeated every time a guess for new 
solution point is required, given the solution at which the 
Jacobian la Invertible. At bifurcation points, results of the 
program DIRECTION are used as a guess for solution in neighborhood 
of the singular point. 


The value Y^ so obtained, is used as guess value to solve the 
following augmented equations by Gauss-Newton ’s method : 

0 = F ( Y ) 

0 = g Y-Y^g - step* (3.18) 

This procedure returns solution as Y^. The tangent vector at 
this point » is found by solving following system of equation. 

0 = E D F( Y’ ) ] <5Y’ (3.19) 

y i 4 


A further check is performed on <5Y’ so as not to take step in 
direction opposite to that taken in previous cycle. Hence, 
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(3.20) 


if. 


<5Y’ . 6Y < 0 

A O 


then the sign of <5Y’ is changed. Y’ is now taken as Y and SY’ as 

A A 0 4 

<5Y^ for the next iteration. The curve length is incremented by step 
size. 


3.8 DYNAMIC STABILITY 


The dynamic stability at each of these points is predicted by 
Routh-Hurwitz criteria. The criteria predict the sign of the real 
part of eigenvalues by merely checking some conditions on the 
coefficients in the characteristic polynomial. For the conditions 
defined by. 


> 0 


a a - 

a a 


a a a 

a a A 


d> ^ 3. ot 

a ft 4 o 


> 0 


a 

- a a 

a o 


> 0 


a a 

A O 


a a 
a a 
a o 


> 0 


(3.21) 


the real part of all the eigen values are negative. 

In the program CONTINUATION , record structure is created to 
store the conditions as functions of coefficients of characteristic 
polynomial. As the necessary data structure for evaluation of 
coefficients are already available, the computations for stability 
can easily be carried out. The Jacobian of the function at that 
point is passed to the boolean function STABLE - It computes the 
coefficients of the characteristic polynomial for the given matrix. 
The conditions coded in DYN_ST are verified for the given values of 
the coefficients. If the conditions are satisfied then the 
function assumes value true. 
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CHAPTER 4 


LOCAL ANALYSIS OF BIFURCATION CHARACTERISTICS 


The formulation and analysis of equation of motion for 
string, done in the previous chapters has been based on the local 
study of the equations about certain critical points, for which 
non-linear nature of the problem is prominent. The governing 
equations have been derived considering the geometeric relationship 
between strain and displacement to be a perturbation from that for 
the linear case. The general motive for application of perturbation 
method to the problem is to facilitate the analysis of the problem 
in the neighborhood of resonance frequency of the linearised 
system, as in such a situation linear theory fails to give 
realistic explanation for the system response. 

Local study of bifurcation characteristics shown by singular 
stationary solutions is an extension of same ideas. Such an 
analysis is imperative for prediction and understanding of the 
general nature of the response. The phenomenon of bifurcation is 
intricately related to the singular points of the map describing 
the general non-linear behavior. 

Bifurcation phenomena occur as a result of interactions 
between solutions predicted by the model. Such an interaction, for 
example, may be an exchange of stability between the solutions. As 
a result of such exchange the response characteristics of the 
system change drastically. In case of strings, there exist 
conditions under which for slight variation in the parameters leads 
to planar motion giving way to non-planar whirling. 
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4-.1 


DYNAMIC STABILITY AT REGULAR POINTS 


Stability o£ the systeoi for a given, value stationary 
solution response can be predicted by the analysis o( the Taylor 
series expansion about that response value. Let, 

dx (4-01) 
dt " 


Then with a change in the coordinates, 

y = X - X 

O 

the differential equation (4.01) may be written as. 


dy _ 

= \y + 


dt “ 


(4.02) 


where is an n x n matrix. This matrix can be reduced to Jordan 
canonical form using a transformation. 


z = Py 

dz _ _ ,, . . 

s + P fix.CP 

H = V * G^(z) 


(4.03) 


where P is taken to be such that det( P ) = 1. 

About regular points the Jacobian matrix A^^ non-zero, 

distinct eigenvalues though they may be complex, hence 


d z. 

dt ^ = «i^i ' 


i = 1, 2,- 


n. 


(4.04) 


In equation (4.04) higher order terms have been negl®^^®*^* 
clear that in such cases stability can be determined *^y computing 
the eigenvalues of the matrix A^^^. If all the eigenval^^® have their 
real parts as negative, then each of the z^ have a form showing 
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0 


exponential decay with time. For such cases the solutions 
(or X = X ), are said to be stable. 

O 

Routh - Hurwitz criteria provide an easier procedure to check 
the signs o£ the real parts o£ eigenvalues o£ a given matrix by 
just cloacking £or certain conditions on coe£ficienta of the 
charact ar-ist ic polynomial ( Chapter 3 ). 

A stable solution becomes unstable when at least one of the 
eigenvalues crosses the imaginary axis in the complex plane 
depicting them, from negative real axis side to positive one. If 
they cnoss for non-zero imaginary value, then it can be argued that 
there ane even number of such eigenvalues and the system now 
becomes capable of showing an oscillatory motion in the respective 
eigen directions. This is a more complex phenomenon than that in 
which t:h.6 eigenvalues are real, and is refered to as Hopf 
bifurea-tion. 

Tliis work is an attempt to analyse the cases in which t 
eigenvalues cross imaginary axis at the origin. For these 
local analysis of the problem about the singular point assumes m 
Importance, as negelecting the occurance of Hopf bifurcation, 
change in stability on a particular branch of solutions which h 
invertible Jacobian, is possible only when a singular point is 
encountered. Thus, branches in the bifurcation diagram on whic 
Jacobian is invertible are capable of showing no bifurcation 
behavior other than Hopf bifurcation. 


4-.2 STABILITY AT A SINGULAR POINT 


Aa, apart from occurance of Hopf bifurcation, all stability 
related phenomena occur only at singular points, and stabili y 
branch of solution with invertible Jacobian, is retained as 
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determined at the singular point from which the branch originated, 
a local stability analysis of each of such branches at singular 
points provides major clues about the general behavior of the 
system. The motive, hence, is to determine the stability at each of 
such curves that emanate from a singular point, locally. 


4.3 ORGANISATION OF SOLUTIONS ABOUT THE SINGULAR POINTS 


The singular points form organising centres for the 
solutions. The conditions under which a singular point exists can 
be described by constraints on the space of parameters. This is 
quite apparent from the fact that for any random value of X the 
solution has invertible Jacobian. The order of a singularity may be 
described by the number of equations that need to be satisfied by 
the parameters for its existence . 

These equations define a smooth manifold in the space of 
parameters. Corresponding to these values of parameters, the Taylor 
expansions of the map are elements of the same subset in the set of 
Taylor aeries of all smooth functions. 


Manifolds so defined may Intersect to give singularities of 
higher order. The singular points can themselves be classified, on 
such considerations. Thus, it can be argued that in a map 
/jj^Cx) v , there exist maximal singularities for some values of 
the parameters, as x is an element of W , a smooth manifold. In 
the case of equation of motion of string. 


/ ( 


X ) = v 


V 


dx 

" at 


(4.05) 


where v = 0 fo stationary solutions. For x = 0 and X = 0, the 

singularity is maximal : it has 4-dimensio. ,nal null space. 
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Organisation of singularities of given map resembles a 
tree-structure in which easing of certain restrictions or variation 
of parameters in a certain direction gives rise to singularity of 
lesser order. This process can be continued upto the singular 
solution of zeroth order, f.e, invertible solutions. The analysis 
about singular point of highest order, provides greater 
understanding about the organisation of solutions and a clearer 
picture of the phenomena of bifurcation that the system might 
show. 


4.4. LOCAL STUDY OF SINGULARITIES 


A general problem faced in the study of non-linear problems 
is the presence of multiple solutions in the neighborhood of 
singular points and ways of obtaining all of them. However, strict 
non-linear nature of the map in the neighborhood of these points 
provides a methodology to study such problems by restricting the 
analysis to broader implications. 


As linear changes in the domain and range do not affect the 
zero set of the Taylor series expansions in terms of bifurcation 
characteristics, equivalence classes may be defined on the set of 
all smooth mappings, under an equivalence relation based on such 
characteristics. The maps which have same bifurcation 
characteristics are members of the same equivalence class. 

For such an analysis, the object of study becomes the 
restriction of the function to a small neighborhood of the singular 

points : 

/ : (EN n 

W S D?” 

IP C R** 

Uher. i3 a s»all opan nei*hborhood o£ the singular poihta 


( 4 . 06 ) 
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The technique proposed for analysis of bifurcation behavior 
relies on the fact that, given an expansion about a singular 
point, a simpler approximate map can be obtained which, in the 
space of all smooth maps, is an element of the same equivalence 
class, defined by the equivalence of bifurcation properties. 

The bifurcation related analysis can then be carried out on 
such an equivalent map. Then, a correspondence can be found between 
the bifurcation behavior and the variation of parameters. 

The set of maps that show similar bifurcation behavior in 
vicinity of a singular point can be Identified with algebraic 
concept of g&r-ms. 


4.4-.1 G&x-ms o/ smooth mappings 

Given a point x « iN , consider the set of all smooth mappings 
from domain •U, which is a neighborhood of x in W, to IP s . On 
this set an equivalence relation, ~ , may be introduced . 

Given two such mappings / : — *■ IP and / : nj — ► IP, / ~ / , 

when there exists a neighborhood of x in W depending on / and / 

A %t 

for which the resrtiction itu coi^'cide. The equivalence 

classes under this relation are called smooth germs of mappings 
£N — > CP at X. The elements of same equivalence class are called 
representatives of same germ. 


Thus, it can be claimed that under the equivalence relation 
defined by type of bifurcation behavior, i .s . , equivalence under 
linear transformations, representatives of the same germ show same 
bifurcation charateristics . 

If and are representatives of the same germ, then 
/^(x) = In view of this fact, notation adopted to denote the 
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germ is, / -• <rD4 , x.> — <rip ,y^, where x and y are called source and 
target of the germ respectively. Given two germs, 

/ : <riN , x> — * CfP , yj> and , 
g : (TIP , yJJ — *. <r® , zJ> , 

operation of composition may be defined as, 

fog: COM, x> — » CO. (4.07) 

This definition does not depend on a particular choice of 
representatives. Under operation of composition there may exist an 
inverse if. 


fog = g o f - 1 


(4.08) 


Inverse mapping theorem applied to germs means that a germ is 
invertible if and only if its differential is invertible. The Rank 
of a germ is defined to be that of its differential or Jacobian. 


4.4.2 Thj& CL — &c^\it'oaX&'njc& of g&r'ms 


Two germs / and / are said to be a-equivanl ent when there 

& ft 

exist invertible germs g and h. for which, if 


then. 



CW^, 


CP^, y^:> 


COM , 

& 

X > — > 
ft 

ft ft 

e ■ 

cm . 

ft 

V - 

cm . x:> 

L A 

K : 

<rip . 

2 

V 7" 

CP , y^J 


f^og-Kof^ 


(4.09) 


The basic objects of study are germs f .-COM, x^ — * CIP, yJ> . 

Hence, this is equivalent to / -• CO? , 0^ — »■ CIR*^, 0^- The set of all 

o 

such germs is termed « 
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« has the structure of a vector space induced from that or 
£R**. Denote to be the set of invertible germs <riR“, 0^ — ► <r[R^ , (1> , 

which is a group under operation of composition. A denotes the 

ry.p ^ 

poduct group A X ^ . The action of the group A on « is 

p rt,p n,p 

given as. 


Ctf , KJ> o / = o / O (4.10) 

o 

Given such an action , a-equivalence of x and y, « •« 

r. ,F 

implies that y = g o x for some g s A and the zero set of each 

encompasses same non-linear behavior. Thus, if x and y are 

a-equival ent , then they are so through a pair of invertible germs 

C h, g ) which is an element of the product group A . Ths 

vs # p 

equivalence classes under such an action are called orbits. Ths 
orbit through X is given by definition as the equivalence class 
which contains x. 


A ox={tfox.-tf«A } (4.11) 

vs.p r,,f> 

The problem of studying the types of bifurcation of the zerc 
set of the mapping ( 4.05 ) is identical to that of studying the 
orbits under this action. 


A-.A-3 Unfolding of g&r^ms 

An r-parameter unfolding of f^ ( r”*, 0) — ► ( , 0 ) is th' 

germ F : ( R^x R*", 0 ) — *■ ( R'x R” 0 ) given by the formula : 

F( u, X ) = ( u , / (u, X)) (4.12) 


4-.4-.4- lief ox-mat ion of get-ms 

An r-parameter deformation of / -• ( R ,0) — > ( R^, 0) is 

germ / ; ( R*" x R”* , 0) — ► ( «"* , 0) with /(0,x) = /^ (x)- 
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THEOREM 4 . 1 


Cll], S&ction B.4.iy 


Let G : <r R * 0 — + <r R**, 0 > be a germ of rank r, the] 

there exists an invertible germ K : G r'', o > v C R*^, 0 > foj 

which, r = <3 o h is an r-parameter unfolding of a germ of rank 0. 


Thus, if the map has y-dimensional null space then it is 
a-equivalent to another map which is (p-y )-paramt6r unfolding of a 
germ of rank zero. 


In Chapter 3, methods have been discussed to restrict the 
equations to the null space .For equation of motion of strings, if 
the null space is y-dimensional , then forms denoting direct 
parameterisations for C 4 - y ) solution variables in the 
compliment of the kernel, can be obtained from equation ( 3.14 ). 

^ ) (4.13) 

o o' 


The Taylor expansion of governing equations about a given singular 
point may be partitioned into kernel and its compliment, and by 
making a change in the variables, 

X - X (x^ , ^•) (4.14) 

0101/9 
o ‘ 

the set of equations are a-equivalent to following map : 

(4.15) 

= d’(y^ . . >-). 

The aim of further theory is to reduce the dimensionality of 
the unfolding of . 0. 0). 



4 -. 4-.5 


Th& a.—&<fui‘i>al 0 rLC& of T-pcor-aitistoT- \a\f o\d.trt£s 


Two r-paraineter unfolding F , F of a germ 

i. S 

f:C ER"', 0^ — *•<!' 03 are said to be a-equivalent when there exist 

r-pararoeter unfoldings I ,2 of the germs at zero of identity 

Wl <1 

maps on 0?”*, IR** respectively for which 

Fj (R^'x R^.O) — ^ (R^x R**, 0) 

F: (R'x R^.O) — CR'xR'*, 0) 

& 

1 : (R'x — *• (R’^x R”, 0) 

w 

I : (R’‘xR'^,0) -♦ (R'xR'’, 0) 

F o J = I o F C4.16) 

4 m <1 ft 



An r-parameter unfolding F of / is a-trivial when it is 
a-equivalent to the constant r-parameter unfolding 

( r' x R"', 0) — ♦ ( R”‘ X R*", 0) of / given by (u, x) — ♦ (u, /(x)). 
/ is called a-stable when every unfolding of / is trivial. 


A-AB 


ContCLC t 


Taking graphs of representatives /^: ( “* ( 

/ : ( JN , X ) — > ( P , y^^), which will be smooth manifolds of 
m X V *and tN^ x through points z^ = (x^ , y^ ) and 

z* = (x , y ), contact equivalence is defined to be the pair 
(h, H) of invertible germs for which one has following commutative 

diagram 


C 


C w X 








where is the inclusion mapping and, n is the 
mapping . Hence , H is given by- 


pro jectior 


H(x,y) = (h(x) , &Cx,y)) 

^Cx,y ) = y 

d St 

Thus, and are contact equivalent or k-equivalent when 
exists a contact equivalence (h, H) for which 

Ho (1,/ ) = Cl,/ ) o h. 

i St 

This is shown schematically in the diagrams given below. 



Fig 2 
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(4.17) 
there 

(4.18) 



Thus, two fierrns (0?", 0)--.(£Rp, 0) lie on the same orbit under 
this action if and only if they are k-equivalent . Thus, the problea 
of classifying germs under k-equivalence la that of classifying the 
orbits under above action, the fc-orbits in 


THEOREM 4.2 C I i } , S&c Itort S. 

Let F, F’: (R X IR ,0) *■ ((R*'x R^,0) be r-parameter unfoldings 

of /p » /jpt of rank 0; If F and F’ are a- equivalent then /j, » » 

are k-equivalent. 


4-.4‘.7 Mor-ph. i srrt. 

Let F , G be r and s parameter deformations of /. A 
morphism from F to <3 is a pair (H, >C) , with K a k-equivalence 
of r parameter deformations , and H a change of parameter, for 
which F is k-equivalent to induced deformation G o H~*' under 

4 . 4.8 V&r-Bctl <d&for-matii>n 

A deformation G of f is k-versal when for any deformation F 

there exist a morphism from F to G . A minimal parameter k-versal 

deformation is called k-universal. 

The assignation ^ induces a mapping from 

a-equivalence classes of germs to k-equivalence classes. Given an 
s-parameter deformation / ; ( R® x R*" , 0) — »• ( R^ , 0) which is 
submersive, the zero set Vy = / ^(0) will be a smooth submanifold 
of R® X r“* of codimension q , £ , a hypersurface defined by 

q-equations in a space of s+m dimensions. Ue define the 
restriction n :C v , 0 .> — *■ <r R* . 0> to be the germ at 0 of 

* * m m 

restriction to of the projection germ Or IR x ik — ER . 
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I. 


THEOREH 4.3 cm 


Se-cttort B. 4. 


Let / ; C tR X R , 0) — ► ( CR*^ , 0) be k-versal 

s-parameter deformation of germ /^ , •• C e”*, 03 — <r e"*, 03 of 

rank 0 ; if are k-equivalent then FI , n are a-equivalent . 


Let 

F : (E*' X E*”, 0> 

► (E'x 


, 0) be an r-parameter 

unfolding 

of a germ 


of 

rank 0, given 

by 

FCu, x) 

= (u, /(u, X)). 

To F, 


we associate a 

germ 

D iCE'x 

E^’x e” , 0) 



. 0), given 

by 


(u, w, x) — *■ -V + /(u, x). Thus, is an (r+q) parameter 

submersive deformation . The Seometric connection 

between F and D is that v is just graph ( / ), and n can be 

P O 1> 

P P 

identified with /. 


THEOREM 4.4 C H 3 » t icn. B. 4. 4J> 

F is a a-stable germ if and only if D is a k-versal 

P 

deformation of / . 

P 

THEOREM 4.5 C[i3 , t ion B. 4. 53 

Let F, <3r CE*'x e”*, 0) ► (E*'xE'^, 0) be stable r-parameter 

unfoldings of the germs and of rank 0; if and /g are 

k-equivalent, then F, G are a-equivalent. 



4..5 OF SMOOTH MAPPINGS 


k 

The jet-space J (n, p) means the real vector space of all 
mappings / •• R — each of whose component is a polynomial of 

degree k in the standard coordinates (x^ , ,x ) in R*' with 

zero constant term. The k-jet of a germ ( J**/ ) corresponds to the 

finite initial segment (degree S k ) of the Taylor series about 
source and target. 

In the set C ( R, R*') given a smooth mapping / — > E** , 

positive numbers « and R, and an integer k ^ 0, a /xindany^rLtctl 

y^igKhc,r^?xc,cd. to / comprises of smooth mappings g .-r” — ^ Rp for 

which , for all x « r" with | x | ^ R , one has 

8 j‘‘/(x) - j’^tfCx) 8 S « (4.19) 

with 8 *=* II * fixed norm on the jet space J** (R*' ,R*' ). 

A subset X c C ( r 7 R**) is called , when given any 

smooth mapping / — ► R** , and any fundamental neighborhood ( v ) 

of / one can find smooth maps a .•r’^ — »■ R^ , in x with « v. 


4.6 FIRST ORDER SINGULARITY SETS 


A singular point of smooth mapping / -* R — R is a point 
X « 04 for which the rank of the differential or Jacobian falls 
below its possible maximal value of min(n, p). 


A natural way of distinguishing one singular point from 
another Is by the actual value taken by the rank of the 
differential. To this end, first order singularity sets are 
defined as follows : 


£['•/ = .| X « W : D / has kernel of dimension i- f 


(4.20) 




“ - { 105 . 921 . 

■cc. No. * '■ ‘ 



is a smooth submanifold of J*^(n,p) of codimension 
i(p-n+i). This is the number of equations or conditions that can be 
imposed on a mappin£ to be an element of . 


4.7 CLASSIFICATION OF SHOCfTH MAPPINGS 

For maps / .•r'^ — R^ first order singularity sets JjV could 
be defined as discussed earlier . If these are smooth manifolds 
then, second order sets can be defined on the restrictions to the 
first order ones, Zf = E* (/ | EV ) . This process can be 
continued. This way higher order TKom. *s singxi’Lojrt s&ts of / may 
be obtained. 


The problem of obtaining these submanifolds £[ in 

k 

the jet-space J (n,p) ,such that for generic / .-R*^ —*■ R**, they 

determine the order singularity sets, was solved by 

S€>ax-db/Tbart C i OQ7^ . 


4 . 7.1 P.&CLI itlg&br-a. 


Real algebra consists of structure of a real vector apace V 
together with a bilinear mapping V x V — > V , written as 
(x, y) — *■ x.y and called algebra product. The real vector space <e 

rt 

of all germs /.-tTR^ — ► <r R, of smooth functions is endowed 
with algebra product induced from that on the reals. 

An ideal is a subset of a ring closed under operations of 
addition and algebra product. The algebra product of an element of 
ideal with any element of ring is an element of ideal. In this 
can be identified with the set of all germs in with 0 as target. 
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4.7.8 


Bocur-dmcm. symbol CUJ, Section B.B> 


Let I be a finitely generated 

» be the generators for I 

system of coordinates in € . Given 

ri 

to be the ideal I + I’ , where I' is 
s X a minors of the Jacobian matrix. 


ideal in the algebra , 

» *nd y , , y be the 

* rs 

an integer s i 1 we define A 
the ideal generated by all 


Given an ideal I c following notation is used : 

A*I = A I 

and A , A , ,A are refered to as successive Jacobian 

extensions of the ideal. As a result of this definition following 
inclusions hold. 


I = A^’l c A*l s= A®I c c A*^! 


( 4 . 21 ) 


I is said to be proper, when I | « . the critical Jacobian 

T\ 

extension of I is the last ideal A*" I in the sequence (4.21) which 
is proper.lt in t'Cifh has a critical Jacobian extension A*'*’’’*'!. In 
this way an ascending sequence A*"*!, a’’**'’’*!, ... , can be 

obtained. I is said to have . ) as Boardman symbol. 


Thus, the Boardman symbol 
defined by that of the ideal 



of a germ / :( IR*',0) — ► ( E**,0) is 
generated by the components of 


TBEOREht 4.6 Clil , S&ctiort B.B.2> 

The Boardman symbol of a germ ( ® ,0) — *• ( is contact 
invariant, i.e., if two such germs are k-equivalent then they have 
same Boardman symbol. 


53 


on the 



should 


The first k integers in Boardman symbol of a germ depend only 
k-jet of /. A necessary and sufficient condition for the set 

, , L k „k 

s= J cn, p) to be non-empty is that following conditions 
be satisfied : 


n ^ V > L > 0 

* ft 

^ n - p (4.22) 

if = n - p then i. =1 = = l 

*■ * ft k 


THEOREM 4.7 C£iJ . Sec tion B. 5. B:> 

If the order singularity set j*’- * « • • • ^ J**Cn, p) is 

non-empty then it is a esTJccth. sxihmani/old. of codimension given by 
the following formula. 


(p -n ) “C\ ■ ^ )^Ci , . • • ,^i ) 


(4.23) 


where ^*(1^,..., ij^) denotes the number of the sequences 

(j ,•...!) of integers which satisfy the following criteria ; 

* k 


j. ^^ft ^ ° 


i 

r r 


(4.24) 


for all s :£ r i k, 
j‘‘(n, p). 


j >0. Hence, the 
are called Bcnaardman. 


singul'arity sets 
t^’uhnuan.i /diie o £ 


THEOREM 4.S <H1 , Section E. S. i > 

Any stable mapping / -• 0?’^ R** is generic ( in the sense of 

Boardman ) . 
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4..8 the procedure of analysis 


Study o£ bifurcation involves the analysis of the solutions 
in the neighborhood of singular points. The solution structure 
about a singular point is the property of the type of singularity 
of the map about that point. Further analysis is based on this 
consideration. 

About a singular point, the map can be modified keeping 
non-linear characteristics invariant. Conditions for such 
approximation are specified by contact invariance and subsequently 
Boardman symbols. 

Satisfying these conditions, an equivalent map can be 
constructed on minimal necceaary information. The bifurcation 
related phenomena can be studied on this map and a correspondence 
can be found between the variation of parameters and the type of 
resulting bifurcation diagram. Such analysis also helps in 
providing good guesses for plotting of exact bifurcation diagram 
about that point. 

For one dimensional case, the maps for which Taylor series 
expansions about a point in variable domain have coefficients of 
all the terms upto secon<|(>rder as zero, show bifurcation 
characteristics similar to those of , 

3 

0 = X + ax + ft (4.25) 

This is a minimal representation of solution structure around 
the singular point. The solutions can be represented in three 
dimensions. Presence of any more parameters do not lead to any 
system behavior which Is not explained by this equation. For 
general r-parameter ( X ) deformation of the original Taylor 
series expansion , same representation holds. The bifurcation 
behavior can still be explained by equation (4.25), execpt that 
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scales on (x, «, ft') are paramet erised by ). Due to 

lesser number of dimensions required to present the bifurcation 
characteristics, the manifold described by the equation may be 
plotted. This provides an aid to predict the type of bifurcation 
that are possible with general variation of (X. , . , . . 

A 1 % 

The set of evolution equations for vibration of the string 
can be represented as a mapping 

/; A X IN — »• W 
W C D?”* 

A = C4.26) 

\y = R"* 

The mapping is expanded about the given singular point . The 
map is considered to be a 5-parameter deformation of / « 

Hence, the initial truncated segment of the Taylor series upto 
order k, J /, can now be taken to be as an element of u ( 4 , 4 ) for 
a fixed value of X. For purpose of representation , however, an 
augmented map is considered which is a B-parameter unfolding of the 
map about the singular point. 


r.-AxIN — lAxV (4.27) 

given by (X, x) — »■ (X, / (X, x)) . The map F is smooth and element 

of * 99 - Introducing a-equivalence on this set theorem ( 4.1 ) can 

be applied. The theorem states that if the germ is of rank p 

9 9 

then there exists an invertible germ h.: C R , 0J> — ^ fR , 0 l> 
such that fS = F o A. is a p-parameter unfolding of germ , g, of rank 
zero . 


The rank of the germ depends on the rank of its Jacobian. As 
5 parameters are already defined, rest of the parameters are 
defined by the rank of the map IN ► V. The map g is the 
restriction of the mapping to the null space ® of IN. R is taken 
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as compliment of «J in IN x A . Ue therefore have. 

<3: CPxQ— ►IPxV 

■SCP. <J) (p. «Cp,q)). (4.28) 

As F and <3 ars a-sq„ivalsnt, theorsm (4.2) helps i„ lookiui for 
another a-equlvale„t p-para«eter unfolding using k-equivalence on 
. Thus the problem of studying 6 -parameter deformation of / is 
transformed to p-parameter deformation of g 

®o ■ 

This part of the operation is oerformarf ir. 

is pertormea in program LYAPUNOV* 

The program ANALYSIS contains routines based on theory discussed in 

the previous sections to represent the map in minimal dimensions. 

The methodology of such reduction is to construct an a-equivalent 

graph in the neighborhood of singular point. 

To e a germ 4 ( ir"‘x e? 0) ro*” is associated given 

by («, w, X) -♦ -w + X). Hence, = 0 defines the graph of 

the equation given by, w = tfCet.x). 

Thus, Dg is an p+ify parameter submersive deformation. In this 
case ni>Q : ( V»q, 0 ) -+ ( , 0 ) can be identified withtf. 

Next step involves an assumption that Dy, is a k-versal 

u 

d 6 f OFinSfli i on » thoFO oxisti srogill dof oFnicftions £of 

which Dg is transverse to all the smooth manifolds of J*'(n, p) 
which define the k-equivalence of g . 

Theorem (4.3) is an observation that if two germs D and D 

G G-' 

are k-veraal p+m parameter deformations of germs g and of 

o 

rank zero and, g and are k-equivalent , then and FId- 

00 G u 

i.e-., G and G’ are a-equivalent. 

At this stage the aim of the operations becomes clear, The 
search is for k-universal such that G and G’ are a-equivalent. 

To find k-universal D^, there is a need to study the k-equivalence 
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class defined by 

The statement of theorem ( 4.6 ) points out that the 

Boardman symbol of a germ is contact-invariant. Hence, to classify 

the germ under k-equivalence , there is need to compute the 

Boardman symbol for and to set up the set of equations by which 

^ k 

the germ g^ is a member of that Boardman submanifold of J Cn»P)- 


4-.9 


BOARDMAN SYMBOL FOR A GIVEN MAP OF RANK ZERO 


The Boardman symbol of a germ corresponds to Thom s 
singularity set. The procedure seeks to construct the set of 
equations by which J**/ is a member of some Boardman submanifold m 
^“(4.4). The map is restriction to the 

q-dimensional null space of the map •-R'* — *• R 

Boardman symbol. Proceeding on similar lines as in Chapter 3, 
Boardman symbol is computed. 


Recall ins the attempt In the previoua chapter to look for 
alngular point or bifurcation point, by aatiafylng condltlona on 

the characteriat ic polynomial such aa equating coefficlenta 
characteristic polynomial to zero, the methodology based o 
Boardman'a theory can be considered to be a recursive use of t e 

same procedure. 


Given a smooth singular map 


with nil-potent Jacobian ( 


the rank = 0 ) , ® 

be satisfied. 


^ , g (x) = V, following equations may 


d ( X ) 


= 0 


C P C X )) = ° 

O q * 


c, P c 3C )) “ ° 


(4.29) 
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Uhere the expressions on the left hand side of the equation set 
(ii ) are the coefficients of the characteristic polynomial of the 
q X q Jacobian matrix of x ). This set of expressions is 
obtained by procedure OET_AUX_EQ . It is to be noted that some of 
these equations have already been satisfied in search of 
bifurcation point . The criteria defining the lc~6quival6nce are 
contructed by sequential transfer of equations from set (n), to 
an another set (ixi), satisfying the condition that the linear 
terms in x. The set (m) is originally empty. 

The methodology in creating set (iii) , is as follows. If 
i ^ and linear terms in it , addition of more 
equations to set (m) is stopped. By a transformation on the 
variable domain a variable ( x^ ) is set up to denote that 
expression by a required linear inverbble transformation and the 
expression is stored in set (m). The set (i) is replaced by the 
set (ix) , with the restriction ( x^= 0 ) and next stage of the 
operations are performed on this set with q-1 equations and the 
process la repeated until all variable directions are so 
determined. The maximum of index t in each stage constitute the 
Boardman symbol. The equations in set (xii) define the Boardman 
submanifold . 


4 ..IO MINIMAL UNFOLDING OF THE GERM OF RANK ZERO 


Operation performed as above can be repeated for the map 
z) = V. The z in this equation are the ones transformed by 
the above procedure. Each of the variables ( z ), represents a 
linear equation in the set (xii). The pararaeterisation of these 
equations by A can again be represented by a corresponding 
translation on . The constant terms which occur in the original 
map C4.26) due to small variations in X , are taken care of by 


translation in V. Thus we have , 


y) = u 

u = u(v, X) (4.30) 

y = y( 2 , X) 

The procedure K_EQ_RED performs the operations mentioned 
above given a germ of rank 0. 

Putting u = 0 and y = 0 in the equation set (in), some 

equation are obtained in terms of X. Parameters are also 

transformed based on these equations , say a = « ( X ) . Dimension 
of « depends on the number of equations in the set at y = 0 and 

u = 0. This in turn depends on co-dimension of the Boardman 
submanifold in the jet-space. These reductions are performed by the 
routine UNFOLDING_X>IR in the program ANALYSIS. 

By making transformations in the domain of Dg it can be 
claimed that Dq,(u, oi, y) is a k-universal deformation , as ot is 
the minimal set of parameters that can describe the versal 
unfolding . 

If « can be represented as to span a subset of tR"*- © , then 
it can be claimed that the transformed map is an a-trivial 
representation in 3*^(4, 4). The affect of small variations of 
parameter for such a representation is that of translation of 
the graph in x V . 


4.i1 CHARACTERISTICS OF A-STABLE GERMS 


If an a-stable representation 
representable in four dimensions or 
variation of parameters does 


is found then graph (/) is 
even less. Thus for these cases 
not change bifurcation 
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characteristics. If the map g(x, X) = v has an a-trivial 

representation , then it is possible to reduce the map to following 
form. 


V’ = X’ ) 

V’ = v’(v, X) (4.31) 

x’ = x’ (x, X) 

where vj (v, X) and x^Cx, X) are invertible. 

Thus, it is feasible to plot the graphC g’ ) in V’x « ’ 
Given a value for X , v’ can be computed using equation C4.31.ii) 
taking V = 0. For that value of v’ the possible solutions ( x’ ) 
can be obtained and since, 

x’(x, X) = X’ (4.32) 

Is invertible, all solutions ( x ) for that value of X can be 
obtained. In this way an approximate bifurcation diagram for a 
given perturbation in the conditions can be plotted with respect to 
any combination of parameters. Having obtained the solutions the 
dynamic stability on each of them may be predicted easily. 
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Chapter 5 


CONCLUSIONS 


The study o£ large amplitude oscillations o£ stretched 
strings provides a representative view o£ the non-linear phenomena 
which occur in actual systems. For the case o£ stretched strings 
such phenomena* e.g.^ whirling and jump, occur near linear 
resonance frequencies. These result from the exchange of stability 
that occurs between the solutions depicting different motion 
characteristics . 

The equations governing the motion of string are formulated 
in Chapter 2 with a view to represent these phenomena on a larger 
time scale. The important assumptions made in the process of the 
formulation are that the wave lengths of the modes are 
comparatively much larger than the amplitudes in the transverse 
directions (2.10), and the velocity of sound in the material of 
the string is much greater than the velocities of the transverse 
modes (2.17). The second assumption reduces the number of the 
differential equations from three to two, equations (2.25) and 
(2.26). 

These equations represent the system to a fair degree of 
accuracy. The governing equations, however, become singular for a 
forcing frequency in the neighborhood of the linear natural 
frequencies. For this reason the course of analysis is to first 
expand the equations about the solutions of linearised system. 
Hence, the spatial displacements are approximated by those of the 
linear modes and the method of multiple scales is applied. This 
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variables in the kernel. A step in these directions and the 
corresponding steps in those spanning the compliment, provides a 
good guess for the program CONTINUATION- 

The program CONTINUATION is used to plot the bifurcation 
diagram for a given parameter. This procedure relies on minimising 
a functional (3.18) to obtain the solutions as they vary with the 
parameter. The solutions so obtained are plotted with respect to a 
given parameter directions. 

These programs were run for a representative set of values. A 
bifurcation diagram was plotted for the given conditions with 
respect to pertubation in the frequency of the forcings from the 
linear resonance frequency, (Fig 3) and (Fig 4) . A bifurcation 

diagram was also plot with respect to amplitude of support movement 
(Fig 5). 

The bifurcation diagrams, plot using such algorithms, do not 
provide much detail about bifurcation charact erstics of 
multi-parameter systems for following reasons : 

(i) All possible invertible solutions arising due to 

variation of a parameter may not be obtained. 

(f£) The continuation scheme is based on minimisation of a 

T 

functional, / / . As the determinant of the Jacobian in the 

neighborhood of singular point is small, variation of 

A { o > 

parameter induces a change of high order ( — ^ — — — - — )in the 

determinant, restricting the corresponding step size in the 
variable domain. 

(fit) These procedures are not capable depicting bifurcation 

diagrams in which singular point is absent. Such bifurcation 
diagrams are important to study the occurance of jump 
phenomena. 
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5.1 


THE AIM OF LOCAL ANALYSIS. 


The program ANALYSIS was coded to aid these procedures in 
ploting of approximate bifurcation diagrams for any set of 
conditions and with respect to any given parameter direction, so as 
to predict the types of bifurcations possible. The technique 
followed, is based on the theories and procedures discussed in the 
Chapter 4. The program, given an unfolding of a germ of rank zero, 
relies on the equivalence conditions defined by the Boardman symbol 
to look for an approximate trivial unfolding. 


The assumption made in the process of such computations is 
that the given mapping is a k-versal unfolding. The pseudo 
parameters can be added to the system to make the unfolding 
k-versal and, on a later stage equated to zero. The program aims to 
map the given unfolding to one of the a-stable representations of 
mappings — »• . The graph of this representation can be plotted. 

The graph so defined corresponds, with invertible 
transformations on the coordinates, i.e. range, domain and 
parameters, to all possible bifurcation diagrams in the 
neighborhood of the singular point. Another characteristic of such 
map is that the bifurcation behavior, with respect to a given 
parameter and for given non-zero values of rest of the parameters, 
does not change with small change in parameters and all solutions 


are invertible. 



Obtaining a-stable representation o£ the mapping under 
a- equival ence relation reduces the problem of analysis to a large 
extent. The classification efforts are aided by the fact that there 
exist only seven types of a-stable representations for mappings 
between 4-dimensional spaces, as listed in S&ctiort S.S.e>. 
Study of such singularities around ( * - 0 ) and ( X = 0) can 
completely unravel all the possible bifurcations of the stationary 
solutions that can occur in the system. 

The possible bifurcation characteristics for a given singular 
point correspond to the Boardman symbol. Some simple bifurcation 
properties possible of generic maps R* — are listed below. 




X 

X. 


Ci = 1, 2, -■ ,4) 


(4.33) 


This symbol denotes a regular point and no bifurcation can 
ever be observed for cases when the map around a point is 
representable in this form ( 4.33 ). 


C i i J> 


If" 


y 

1 . 



X 


2 

X 

4 


(i 


1, 2, 3) 

(4.34) 


The maps a-equivalent to above map can show bifurcation 
characteristics in which two solutions emerge from the point. 
Bifurcations represented by a pair of intersecting lines are 
represented by forms of this type. For the germs of this class the 
solution may disappear for some value of parameters. 





y. 


X 


a 

X 

<4 


+ X X 

A ^ 


(i 


1, 2, 3) 


(4.35) 


The singular points which have the form (4.35) associated 
with them are capable of displaying pitch-fork bifurcations and the 
bifurcation characteristics arising due to breaking of the 
pitch-fork. The singular point obtained for (Fig 3) is one of thie 
type. It is to be noted that for such cases at least one solution 
is always present as against previous case when for some value of 
parameters solution could have disappeared. These cases ere 
important for general stability analysis as phenomena of jump oan 
first be observed at only such points. 


<r 




y. = X. 


(i = 1, 2, 3) 


= X 


X X 
4 


2 

X X 
2 -A 


(4.36) 




j^,4,4^4,0 


y. = x^ 


(i = 1. 2, 3) 


= X 


X X 

4 4 


2 

X X 
2 -A 


2 


(4.37) 


These cases represent more complex solution organisation 
around singular points of the respective types, from which upto 
four ( 4.36 ) ov five C 4.37 ) solutions may emerge on variation of 

certain paramet ers . 





y 5=x. (i=l,2) 

i 1. 

y = X X 

a -A 

y = X* + X* + X X^ + X X 
fl .A 42 2-A 


(4.38) 
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Sin-gulctr- Poirtt 1.4966 0.310 0.0 0.0 0.01 3.0 1.051 5.0 116.68 



TSOf’i.BinOBO 



y = X 


(i = 1, 2) 






y. = 


X X 

a •« 

2 2 

X - X 

a * 


X X 

£ a 


X X 
8 * 


(4.39) 


Both these forms represent more diverse bifurcation phenomena 
than the previous cases. Form (4.38 ) represents the solutions 
arising from intersection of a hyperbola and a circle in a plane 
defined by x^ and x^ . The expansions about singular points which 
are a-equivalent to form ( 4.39 ) show bifurcations arising from 
intersection of two hyperbolas. 

The program ANALYSIS was coded for the purpose of obtaining 
such representations for given expansions and analysing them. 
However satisfactory runs were not possible due'bo irregularity and 
unavailability of computational facilities and time. 
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